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Abstract 

An analysis of the existing data on the atmospheric neutrino anomaly is pre- 
sented, focused on the statistical significance that can be attributed to its 
experimental evidence. Our approach is alternative to the usual analyses in 
terms of the fi/e ratio of event rates. In fact, we perform a comparison be- 
tween data and expectations, by separating the information on e-like and /x-like 
events, with a careful estimate of the different errors and of their correlation 
effects. The results are shown both numerically and graphically, and disclose 
interesting aspects of the atmospheric neutrino anomaly, that the use of the 
fi/e ratio would partially hide, both in the sub-GeV and in the multi-GeV 
energy range. 

1 Introduction 

"Atmospheric neutrino anomaly" is usually referred to as the unexpected difference be- 
tween measured and predicted muon/electron flavor composition of the atmospheric neu- 
trino flux. Claimed as possible evidence of new physics beyond the Standard Model of 
electroweak interactions, it is generally interpreted in terms of neutrino oscillations. 

First pointed out by the Kamiokande collaboration [1], the evidence for an anomaly 
in the sub-GeV energy range ((E v ) ^ 1 GeV) has been reinforced in further exposures of 
the same Kamiokande detector [2, 3], and has also been confirmed by a similar (water- 
Cherenkov) underground experiment, 1MB [4, 5]. More recently, the Kamiokande collab- 
oration has reported that the flavor composition of a higher energy event sample (the 



so-called multi-GeV events) is also anomalous [3]. On the other hand, two of the iron- 
calorimeter experiments, Frejus [6, 7] and NUSEX [8], did not find results in conflict with 
the expectations. The third, Soudan 2 [9], possibly does, although its data analysis is still 
preliminary. 

Actually, a comparison of the experimental data with the expectations requires a 
reliable and precise calculation of the (anti)neutrino fluxes and their flavor composition. 
Conversely, as is well known, we observe a large spread among the different independent 
atmospheric v flux calculations, hereafter referred to as: BGS (Barr, Gaisser and Stanev) 
[10], P (Perkins) [11], HKHM (Honda et al.) [12], LK (Lee and Koh) [13], BN (Bugaev 
and Naumov) [14], KM (Kawasaki and Mizuta) [15] 1 . This spread reflects essentially 
the large uncertainty in the overall normalization (of order 20-30%). The uncertainty, 
however, is reduced to a value as low as 5% when the \ije flavor ratio is considered [17]. 
This cancellation of errors is the main motivation for the commonly used double ratio 

R»/e = (fJ-/e) D ata/(fJ-/e) M C (1) 

(see, e.g., the systematic review [18]) where /x, e represent, respectively, the number of 
/x-like and e-like events in a given detector, as observed (Data) or simulated through a 
Monte Carlo (MC) numerical experiment. 

The (often forgotten) drawback of the above double ratio is that, by construction, its 
error distribution is non-gaussian. A typical "anomalous" result, such as R^/ e = 0.5 ±0.1, 
represents thus an uncomplete summary of the data that could mislead the naive reader, 
if not supplemented by the true error distribution. Such distribution should be given 
explicitly by the experimental collaborations, since it involves the knowledge of the errors 
before the ratio is taken. The use of two slightly asymmetric errors on R^/e as in Ref. [3], 
is only of limited help, and certainly is not sufficient to recover, for people other than the 
experimentalists themselves, the full information needed for correct statistical tests and 
phenomenological analyses. 

Our point of view is that it is possible (and perhaps easier) to use exclusively gaussian 
distributions, and also to exploit fully the experimental and theoretical information, by 
separating the e and \i flavor data. The modest price to pay is that the correlation of the 
different variables must be accounted for. More precisely, we will use and compare only 
those variables whose errors can be assumed to be normal (at least in the absence of any 
a priori contrary reason), that is: fi M c, fJ-Data, e M c and e Data . 

In the following sections this approach will be presented and developed systematically, 
in order to assess the significance of the atmospheric neutrino anomaly, as revealed (or 
not) in the various experiments performed so far. More precisely, in Section 2 we discuss 
in detail the experimental and simulated data, with a careful analysis of their errors and 
correlation effects. Section 3 is devoted to a comparative discussion of the results, with a 

1 All these neutrino flux calculations include the muon polarization effect, first pointed out by Volkova 
[16]. 
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specific example of a possible interpretation in terms of neutrino oscillations. In Section 4 
we consider the very interesting analysis of the multi-GeV events, binned according to 
their direction. Finally, in Section 5 we draw our conclusions 2 . 

2 Analysis of (/i, e)r>ata an d (/i, e)MG uncertainties 

In this Section the gaussian errors affecting the doublet of variables (/x, e)u ata and (/x, e)MC 
are estimated together with their correlations. We discuss first the largest, common, 
source of uncertainty in (/x, e)MC, induced by the spread in the theoretical calculation of 
the neutrino fluxes. Then, we estimate all the other relevant error sources for each specific 
experiment. 

As far as neutrino fluxes are concerned, let us consider the predicted (/x, e) rates for the 
Kamiokande detector in the sub-GeV range, as calculated by using the BGS input fluxes 
(from this point forward, we conventionally rescale all rates to the central value of the 
BGS rates). Accordingly, the theoretical predictions, treated as a statistical population, 
can be conveniently represented by a bivariate gaussian distribution in (/x/ix^GSj e/e^Gs) 
with the center at (1, 1), errors s e and s M as large as 30% and correlation p = 0.986. The 
standard deviation ellipse 3 corresponding to such distribution is shown in Fig. 1. In the 
same figure we also show the alternative predictions (only the central values) coming from 
the other different input fluxes (P, HKHM, LK, BN, KM). In particular, the spread in 
the KM predictions (dots) is reminiscent of the flux variations obtained by the authors 
by varying a few input parameters [15]. Figures for other detectors (not shown) would be 
similar. In Fig. 1 we also draw some isolines of the double ratio (xx/e)/(xx/e)£GS (dashed), 
that will prove useful in the following. 

Let us briefly comment some of our previous choices: 

1) We choose to center the distribution on the BGS predictions because the BGS neutrino 

spectra are very well documented and detailed, and have also been used by all the 
experimental collaborations in at least one simulation, so that they are appropriate 
for a global comparison. 

2) A la uncertainty as large as 30% accounts conservatively for those flux calculations 

having the smallest normalization (LK and BN), and, at the same time, reduces 
much of the model-dependence implicit in the previous choice of BGS as reference 
fluxes. The correlation value p = 0.986 guarantees a residual s M / e = 5% theoretical 
error on the /x/e ratio. 

2 This work does not cover the phenomenology of neutrino-induced upward-going muons, for which 
there are already extensive analyses [19, 20] as well as critical insights [21, 22]. 

3 We use the notation of Ref. [23]. The standard deviation ellipse corresponds to the 39% C.L. in two 
variables; its projections onto the axes correspond to ±lcr errors on e, /x. The ellipse containing 68% 
(90%) probability would be similar, but rescaled linearly by a factor 1.51 (2.15). 
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3) In principle, the flux uncertainty could be substantially reduced [17] by using as an 
additional constraint the data on negative muon fluxes, as those reported in [24]. 
This constraint is used explicitly in Ref. [11] (which is in good agreement with BGS) 
and in a very recent and detailed flux calculation [25]. In particular, in the latter 
work the value s M / e = 5% is confirmed, and the estimate = s e = 20% is defended. 
We will thus supplement our "default" case (sfi ux = 30%, used in all figures) with 
the less pessimistic error estimate Sfi ux = 20%. In this case, the corresponding 
correlation is p = 0.969. The purely hypothetical case = s e = 10% will be also 
discussed. 

A further comment on Fig. 1 is in order. A large cancellation of the theoretical flux 
uncertainties down to 5% can be considered established only for the ratio of the total 
rates, but is not guaranteed in small subsamples of the MC events, as, e.g., the five bins 
of the angular distribution of multi-GeV events (Ref. [3], Fig. 4), that will be studied 
in Sec. 4. Indeed, as far as the theoretical angular flux distributions shown in Ref. [15] 
(HKHM) and Ref. [26] (KM) can be compared, the bin-by-bin differences in v^j ' v e appear 
to be larger than 5%. The same trend, i.e. larger uncertainties in each bin than in the 
integrated rate, characterizes the fluxes at the higher energies relevant for upward- 
going muon production [19, 27]. Lacking a detailed comparison of the different predicted 
neutrino angular distributions, it seems reasonable to assume an uncertainty of 10% on 
the fi/e ratio in each of the above five bin, this choice being more conservative than that 
performed in the Kamiokande analysis [3]. 

The above theoretical uncertainties are shared by all the experiments. Now we have 
to consider the remaining detector-dependent errors, which affect both (/x, e)u ata and 
(Mj e)*fc in each different experiment. Thus, we list in Table 1 the absolute (published) 
values of (/x, e)u ata and (/x, e)MC for the various atmospheric neutrino experiments. For 
the reasons given before, all the simulated numbers correspond to BGS input fluxes. Also 
shown in Table 1 are the total and simulated exposures for each detector. Concerning 
the Kamiokande experiment, the reported data refer to both the sub-GeV and the multi- 
GeV energy range, the latter including fully contained (FC) and partially contained (PC) 
events [3]. The FC and PC samples are characterized by (E v ) ~ 3 GeV and (E v ) ~ 9 GeV, 
respectively. In the multi-GeV range, the quoted number of MC events corresponds, more 
precisely, to BGS fluxes supplemented with Volkova fluxes [28] for E v ^ few GeV ("flux 
B" of Ref. [3]). The additional information provided by the zenith-angle distribution of 
multi-GeV data will be examined separately in Sec. 4. For the 1MB detector [4, 5], only 
the contained event sample is considered. The data for Frejus in Table 1 include the fully 
contained (FC) events, and the total sample (ALL), as reported in [6, 7]. Concerning 
the NUSEX experiment, the predictions apparently refer [8] to the BGS fluxes without 
muon polarization [29]. In view of the very large errors in this experiment, we have not 
attempted any correction. Our source for the NUSEX simulated exposure is Ref. [30]. In 
the last row of Table 1 we consider the preliminary Soudan 2 data [9]: in this case, the 
number of observed events is not an integer, due to a "shield inefficiency correction" [9]. 
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We now present our analysis of the uncertainties. Although only published data are 
used (unless otherwise noticed), our method differs from the usual approaches, in that we 
are interested in the errors affecting the separate flavors p and e, and not only the flavor 
ratio Rft/e . This requires that the error correlations p^ e are to be taken into account, as 
done previously for the flux uncertainties. The reader is referred to Table 2, where we 
have collected the actual values of the individual sources of errors and their combinations. 

In the discussion of Table 2, let us first consider the data errors. Of course, all data 
samples are affected by statistical errors (with p^ e = 0), obtainable from the values of 
l^Data and euata of Table 1. All the experiments also consider the possibility of flavor 
misidentification (mis-id), for which we have taken the published values, with full anti- 
correlation: p^ e = — 1. Kamiokande lists several additional sources of data errors: multi- 
ring event separation, vertex fit, absolute energy calibration and non-neutrino background. 
For sub-GeV data, they add up to s M / e ~ 2.7%. In the absence of any other information, 
in Table 2 their correlation is disregarded and they are assumed as equally shared between 
p and e: s e = = 2%. Analogously, for multi-GeV data it is s M / e ~ 5.4% [3], so that 
s e = = 3.8%. Concerning the Frejus experiment, the trigger efficiency uncertainties 
[6] should also be considered. They can be disregarded for p-like events, but are sizeable 
for e-like events: s e = 10%. No other sources of data errors are quantitatively discussed 
in the NUSEX paper [8] and in the Soudan 2 report [9]. Finally, the "Total" errors in 
Table 2 are obtained by summing in quadrature all the 2x2 error matrices associated to 
the values of s M , s e and p^ e considered so far, for each experiment separately. 

Let us now consider the MC errors, i.e. those affecting the simulation of event pro- 
duction in each detector. The largest contribution, provided by the neutrino flux uncer- 
tainties, has already been discussed at length, and is not reported in Table 2. For all 
the simulated samples, the statistical errors must be taken into account, according to a 
binomial distribution of p and e events. The relevant input can be taken from Table 1. An 
explicit check of our estimate is possible for Kamiokande, the quoted statistical MC errors 
being s M / e = 3.6% [26] and s M / e = 6% [3] for the sub-GeV and FC+PC multi-GeV cases, 
respectively. Our estimates in Table 2 imply that s^/ e = 3.4% and s^/ e = 5.6% respectively 
(the agreement would be even better by using HKHM instead of BGS fluxes). Concerning 
the MC errors related to the neutrino interaction in the detectors, cross section and nu- 
clear model uncertainties have been treated differently by the various collaborations. For 
Kamiokande sub-GeV, charged current (CC) cross section errors are estimated to amount 
to ~ 10% for each flavor [31], reduced to 3% in the p/e ratio [26], p^ e being determined 
by this cancellation. Neutral current (NC) cross section uncertainties affect mainly e^c 
through 7T° contamination in the Kamiokande sample: s e ~ 2%, s M ~ 0. These numbers 
are slightly different for the multi-GeV case: s M / e (CC) = 2% and s e (NC) = 3% [3]. The 
1MB collaboration estimates nuclear and cross section uncertainties more conservatively 
than Kamiokande, although it is difficult to extract definite error values from the pub- 
lished papers [4, 5]. The single largest effect in the 1MB simulation is induced by varying 
by ±20% the axial mass parameter, leading to s M ~ s e ~ 20% and s M / e ~ 10% [32]. Fermi 
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gas model uncertainties are estimated to increase s^/ e up to ~ 14%. The final values of 
"nucl.+cross" in Table 2 for 1MB correspond to the choice (s M , s e , s M / e ) = (20, 20, 14). 
For the Frejus experiment, nuclear and cross-section uncertainties amount to s M / e = 6% 
[7]. Values of s M and s e are not published, however taking s M ~ s e ~ 10% (similar to the 
Kamiokande case) is not unreasonable, and leads to the values in Table 2. Concerning the 
NUSEX and Soudan 2 experiments, in the absence of detailed published information we 
have assumed nuclear uncertainties as large as for Frejus, a choice which hardly leads to 
overestimate the errors. Finally, the "Total" MC errors reported in Table 2 are obtained 
by adding in quadrature the previous errors and the flux uncertainties (only the case 
Sfiux = 30% is explicitly reported). 



3 Synopsis of the Results 

We have estimated, as carefully as we could, the errors affecting (/x, e)u ata and (/x, e)MC 
for each experiment. There are no a priori reasons against the assignment of gaussian 
distributions to these errors, except perhaps for the fluctuations in the (small) number 
of events observed by NUSEX, which should be described more properly by a Poisson 
distribution. However, data and simulations do agree in NUSEX, so that the tail of 
the statistical error distribution is not probed, and the Poisson distribution can be well 
approximated by a gaussian within ±lo\ 

Thus, an unbiased % 2 can be defined for each experiment: 

X 2 = E ^■(<T 2 Data + <T 2 MC )~ a 1 p -Ap , (2) 

a, (i=n, e 

where a 2 are the 2x2 total (squared) error matrices, and A is the vector of the Data— MC 
differences: A = (fiData ~ fJ-MC, e Data - e M c)- 

The list of % 2 values and corresponding C.L. (d.o.f. = 2) for the experiments examined 
in the previous Section is given in Table 3. The second and third column refer to our 
"default" case (sfi ux = 30%), the next two columns to the estimate Sfi ux = 20% reported 
in [25] and the last two columns to the purely hypothetical case Sfi ux = 10%. In any given 
column, the "hierarchy" of experimental evidences for the atmospheric neutrino anomaly 
is rather evident and does not need any comment. The comparison between the columns 
at Sfi ux = 30% and 20% shows that the results are remarkably stable with respect to the 
normalization error, especially for the experiments that do show an anomaly, the reason 
being that, in this case, the greatest contribution to the % 2 is related to the uncertainties 
orthogonal to Sfi ux . A further hypothetical reduction down to 10% breaks this stability, 
increasing a few % 2 values much more than others; in this situation, however, also the 
specific choice of a "reference" flux (BGS in our case) starts becoming crucial for the 
results. 
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The above results can be better understood by means of Fig. 2, where the standard 
error ellipses (39% C.L.) of data (white) and MC (gray) are displayed for each experi- 
ment, together with the corresponding iso-i? M / e lines (dashed). All MC ellipse include a 
30% flux error. The 1MB and Soudan 2 MC ellipses are the largest ones, as a result of 
more conservative nuclear error estimates (1MB) or small simulated statistics (Soudan 2). 
The narrowness of the NUSEX Monte Carlo ellipse only reflects our ignorance of the 
corresponding systematics. 

Before commenting on the single experiments shown in Fig. 2, it should be noted that 
any change in the overall normalization of the fluxes (e.g., any choice of reference fluxes 
different from BGS) would only have the effect of shifting the MC gray ellipses up or down 
along their major axis. Moreover, such shifts are bound to be approximately equal for all 
those experiments which are sensitive to the same range of the neutrino energy spectrum. 
In particular, this property holds for NUSEX, Frejus (FC), Kamiokande sub-GeV, 1MB, 
and possibly Soudan 2: all of them, in fact, observe contained events and are sensitive to 
0.2 <E V <2 GeV. 

Concerning the iron-calorimeter experiments (Fig 2a,b,c), the close agreement be- 
tween the NUSEX and Frejus results on \i and e separately is remarkable, being more 
informative than simply the agreement of the R^/ e values (dashed lines). Both exper- 
iments seem to favor fluxes with low normalization. The preliminary Soudan 2 data 
are slightly far from the expectations. Fig. 2c provides the additional information that 
Soudan 2 data favor fluxes with "central normalization". Assuming a theoretical error 
lower than 30% in the MC ellipses, the different indications provided by the three iron 
detectors would be correspondingly exacerbated. 

Concerning water- Cherenkov experiments (Fig. 2d,e,f), it is the very good agreement 
between the data of the two high-statistics experiments 1MB and Kamiokande sub-GeV 
(white ellipses), and their common disagreement with the MC simulations (gray ellipses), 
which provides the well-known evidence for an anomaly in the sub-GeV range. It is 
interesting to observe, however, that the standard deviation error ellipses of the 1MB, 
Kamiokande sub-GeV, Frejus FC, NUSEX and Soudan 2 data are mutually compatible. 
The relative position of these five data ellipses would not be spoiled by choosing a reference 
flux different from BGS: as said, that would simply correspond to shift the five MC ellipses 
by one and the same amount. Thus, this additional degree of freedom cannot bring to a 
closer agreement the indications coming from Frejus and NUSEX on the one hand, and 
Kamiokande sub-GeV and 1MB on the other hand. 

Multi-GeV data (Fig. 2f) show an additional feature. If we limit our attention only 
to the ratio -R M / e , the results from the sub- and multi-GeV samples of Kamiokande agree 
rather well (slanted lines of Figs. 2e,f). But this is not the case for the separate \i and e 
flavours, as it clearly emerges from the different position of the data ellipses. This does 
not imply an inconsistency between the two data sets: simply, it seems to indicate that 
the ratios \ij \ibgs and e/e^GS increase with energy, while keeping the double ratio R^/ e 
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approximately constant. We note that part of this effect could be explained, for instance, 
by assuming a corresponding hypothetical decrease in the slope of the theoretical neutrino 
energy spectra. With regard to this, future re-calculations of atmospheric fluxes can 
usefully address the problem of estimating the allowed range of such shape variations. 

The more abundant information coming from the comparison of sub- and multi-GeV 
data (not even including the directional information of the multi-GeV data) can make the 
usual attempts to explain the atmospheric neutrino anomaly less feasible. As an example, 
we show in Fig. 3 what happens by assuming pure — > v T oscillations with a high value 
of Am 2 , so that P^v^ — > v T ) — \ sin 2 29, independently on the energy (the e flavor is 
not affected at all). All error ellipses are considered in the same plot, whereas, in order 
to avoid confusion, only a representative MC ellipse is shown, both in the standard case 
(dashed) and assuming oscillations (gray). The choice sin 2 29 = 0.66 is seen to bring the 
theoretical predictions in agreement with the data coming from the different experiment. 
However, the sub- and multi-GeV Kamiokande data clearly favor different normalization 
of the MC rates. With theoretical errors much smaller than 30%, — > v T oscillations 
with high Am 2 would thus be unable to fit both sub- and multi-GeV data at the same 
time: this could be a tentative indication in favor of additional effects, able to renormalize 
the rates in an energy-dependent way. 

Apart from the above example, and another sketchy analysis of <-> v T mixing in 
the next Section, we do not further investigate in this paper neutrino oscillations, as an 
exhaustive analysis of them requires, in our opinion, a separate work [33]. 

4 Analysing Multi-GeV Binned Data 

So far we have analyzed only total rates, and not the full information contained in con- 
venient subsamples, as the histograms reporting energy and/or angle distributions. It 
is difficult to perform accurate "binned analyses," for at least two reasons. On the one 
hand, bin-by-bin experimental systematics are usually not published, and cannot even be 
guessed. One the other hand, theoretical uncertainties are less under control, since their 
dependence on the neutrino energy and/or direction is largely unknown. 

Nevertheless, a few interesting insights can be gained by looking at a specific example: 
the angular distribution of multi-GeV data in Kamiokande [3]. This is a particularly 
important case, since the evidence for an anomalous angular distribution would strengthen 
the neutrino oscillation hypothesis. We use the information contained in Figs. 3a and 3b 
of Ref. [3], where the data were divided into five bins equally spaced in the cosine of the 
zenith angle. The first (fifth) bin correspond to upward (downward) neutrinos. In each 
bin, we compare (//, e)u ata with (//, e)MC as it has been done before for the total rates, 
the only exception being that now we allow (.s M / e )/z U x = 10%, as discussed in Sec. 2. 

The results are shown in the upper part of Fig. 4, reported as the "no-oscillation case," 
with the uncertainties drawn as gray (MC) and white (Data) ellipses. Data errors are 
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assumed gaussian, although a Poisson distribution would be more appropriate to describe 
small data samples. This choice, however, would imply making the variables (/x, e)u ata 
discrete, thus preventing any simple graphical representation. Moreover, it can easily 
be checked that for the integer values of (/x, e)u ata of interest in Fig. 4, the difference 
between poissonian and gaussian ±lcr intervals is small, and in any case irrelevant for 
the following discussion. Coming back then to this "no-oscillation case," we see that the 
agreement between data and MC simulation, quite bad in the 1st bin, increases in the 
next four bins. 

A better fit can be obtained, however, by assuming a neutrino oscillation scenario. 
Again, the simple case of pure — > v T oscillations is considered. In particular, the 
medium strip of Fig. 4 displays the results obtained for maximal mixing (sin 2 29 = 1) 
and large Am 2 (Am 2 > 10 _1 eV 2 ), so that \imc —> \^mc- In this case, the agreement 
between data and MC simulation is lost in the 5th bin, but improves in the first 4 bins. 

An even better fit can be obtained by lowering Am 2 . In particular, the lower part of 
Fig. 4 shows the analysis for the best fit values of the mass/mixing parameters, as taken 
from Ref. [3]. Now the agreement between data and MC simulation is improved with 
respect to both the previous cases. It is instructive to observe that the simple comparison 
of fi/e ratios, corresponding to the slopes of the slanted lines, would have hidden the 
persistent discrepancy between data and MC simulation in the first bin. 

We do not attach any definite C.L. to the three scenarios shown in Fig. 4, because of 
the aforementioned ignorance of potentially important bin-by-bin correlations, and leave 
the reader to judge the significance of the angular anomaly and its explanation in terms 
of oscillations. As far as our opinion is concerned, we quite reasonably conclude that the 
evidence for the third (oscillation) scenario, although significant, is not really striking. In 
particular, it should be noted that the 1st and the 5th bin, which play an important role 
in the fit, contain the smallest number of events, the angular \i and e distributions being 
both peaked at the central bin [3]. 

It must be noted that the information used by the Kamiokande collaboration for statis- 
tical tests is larger than the angular distribution alone, and includes a 8 X 5 (energy X angle) 
histogram of MC and Data events. Unfortunately, only the energy and angle projections 
are published [3]. We hope that in future publications this additional, and potentially 
very important, spectral information will be fully reported. 

5 Conclusions 

The atmospheric neutrino data coming from a large part of the experiments performed so 
far show an interesting pattern of deviations with respect to the theoretical predictions, 
that is usually summarized in the double ratio R^/ e . This requires, however, an analysis 
of the non-gaussian distribution of the ratio uncertainties. 
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An alternative approach has been proposed in this paper, which allows to use only 
gaussian distributions, provided that the information on the the \i and e flavors are 
separated, and their correlation effects are properly taken into account. 

Accordingly, we have performed for each experiment a careful analysis of the corre- 
lated uncertainties affecting the e-like and mu-like observed and simulated event rates. 
This enabled us to assess quantitatively the statistical significance of the atmospheric 
neutrino anomaly in each experiment, and to study the sensitivity of the results to the 
flux uncertainties. 

The anomalous angular distribution of the Kamiokande multi-GeV data has been also 
analyzed, separately, with the same methodology. It is shown graphically that neutrino 
oscillations can bring the expectations in closer agreement with the data, although more 
data are needed to derive compelling indications. 

In conclusion, we have shown in detail an unconventional way of looking at the atmo- 
spheric neutrino anomaly, in which the slight complication of taking into account non-zero 
error correlations is more than compensated by the benefit of having clear graphical and 
numerical results. 
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Figure Captions 



Fig. 1: Comparison between the different theoretical predictions of the rates of /x-like 
and e-like events (in the sub-GeV range) for the Kamiokande detector, normalized 
to the predictions of the reference neutrino fluxes of Ref. [10] (BGS fluxes). The 
standard deviation ellipse (39% C.L.) assumed to represent the maximum (30%) 
theoretical dispersion is also shown. See the text for details. 

Fig. 2: Comparison between Monte Carlo simulations (gray ellipses) and experimental 
data (white ellipses) for the following experiments: NUSEX, Frejus (ALL events and 
only fully contained, FC, events), Soudan 2 (preliminary), 1MB, Kamiokande sub- 
GeV, Kamiokande multi-GeV (fully contained, FC, and partially contained, PC). 
Iso-lines of the double ratio R^/e are also shown (dashed). 

Fig. 3: Simultaneous comparison of the data with the theoretical (MC) predictions, 
including an example of neutrino oscillation effects. Dashed ellipse: MC, without 
neutrino oscillations. Gray ellipse: MC, with — > v T oscillations at large Am 2 and 
sin 2 26 = 0.66. 

Fig. 4: Comparison of multi-GeV Kamiokande data (white ellipses) and simulations 
(gray ellipses) in each of the five zenith-angle bins. Upper strip: no neutrino oscil- 
lations. Medium strip: — > v T oscillations with maximal mixing and large Am 2 . 
Lower strip: — > v T oscillations with best fit mass/mixing parameters. 
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Table 1: Number of /z-like and e-like events, as observed in each experiment (Data) or 
simulated with the reference fluxes of Ref. [10] (BGS). Event numbers refer to the total 
exposure of the detector (the exposure used in the MC simulation is also given). 







Total 


Simulated 










Experiment 


Refs. 


exposure 


exposure 


/x-like events 


e-like events 






(kt-yr) 


(kt-yr) 


Data 


BGS 


Data 


BGS 


Kamiokande 
sub-GeV 


[3] 


7.7 


43 


234 


396.0 


248 


257.2 


Kamiokande 
multi-GeV (FC) 


[3] 


8.2 


51 


31 


40.4 


98 


70.8 


Kamiokande 
multi-GeV (FC+PC) 


[3] 


8.2+6.0 


51+40 


135 


165.8 


98 


70.8 


1MB 


[4] ,[5] 


7.7 


33 


182 


344.5 


325 


339.4 


Frejus (FC) 


[6] ,[7] 


1.56 


10 


66 


90.0 


56 


66.8 


Frejus (ALL) 


[6] ,[7] 


1.56 


10 


108 


125.8 


57 


70.6 


NUSEX 


[8] 


0.74 


~15 


32 


36.8 


18 


20.5 


Soudan 2 
(preliminary) 


[9] 


1.01 


3.72 


33.5 


42.1 


35.3 


28.7 
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Table 2: Individual and combined errors s M and s e (%), and their correlation p^ e , for 
both observed (Data) and simulated (MC) samples of /x-like and e-like events. Total 
MC errors implicitly include neutrino flux uncertainties, characterized by: (s M , s e , p^ e ) = 
(30,30,0.986). 



Experiment 


Data errors 
Source 5 M s e p^ 


MC errors 
Source s e p^ 


Kamiokande 
sub-GeV 


statistics 6.5 6.4 0.000 
mis-id. 2.0 2.0 -1.000 

various 2.0 2.0 0.000 


statistics 1.3 2.1 -1.000 
CC cross sec. 10.0 10.0 0.955 
NC cross sec. 0.0 2.0 0.000 


Total 7.1 7.0 -0.081 


Total 31.7 31.8 0.975 


Kamiokande 
multi-GeV 
(FC, FC+PC) 


stat.(FC) 18.0 10.1 0.000 
stat. (FC+PC) 8.6 10.1 0.000 

mis-id. 2.0 2.0 -1.000 

various 3.8 3.8 0.000 


stat.(FC) 5.0 2.9 -1.000 
stat. (FC+PC) 1.6 4.0 -1.000 
CC cross sec. 10.0 10.0 0.980 
NC cross sec. 0.0 3.0 0.000 


Tot.(FC) 18.5 11.0 -0.020 
Tot.(FC+PC) 9.6 11.0 -0.038 


Tot.(FC) 32.0 31.9 0.951 
Tot.(FC+PC) 31.7 32.0 0.966 


1MB 


statistics 7.4 5.5 0.000 

mis-id. 5.0 5.0 -1.000 


statistics 1.8 1.9 -1.000 
nucl.+cross 20.0 20.0 0.755 


m , i on n c n on a 
Total o.y (.0 —0.3(4 


Total 00. 1 OD. 1 O.yiO 


Frejus 
(FC, ALL) 


stat.(FC) 12.3 13.4 0.000 
stat.(ALL) 9.6 13.2 0.000 
mis-id. 2.0 2.0 -1.000 

trigger eff. 0.0 10.0 0.000 


stat.(FC) 2.7 3.7 -1.000 
stat.(ALL) 2.1 3.8 -1.000 
nucl.+cross 10.0 10.0 0.820 


Tot.(FC) 12.5 16.8 -0.019 
Tot. (ALL) 9.8 16.7 -0.024 


Tot.(FC) 31.7 31.8 0.950 
Tot. (ALL) 31.7 31.9 0.953 


NUSEX 


statistics 17.7 23.6 0.000 
mis-id. 5.0 5.0 -1.000 


statistics 2.2 4.0 -1.000 
nucl.+cross 10.0 10.0 0.820 


Total 18.4 24.1 -0.057 


Total 31.7 31.8 0.951 


Soudan 2 
(prelimin.) 


statistics 17.3 16.8 0.000 
mis-id. 5.0 5.0 -1.000 


statistics 5.0 7.5 -1.000 
nucl.+cross 10.0 10.0 0.820 


Total 18.0 17.6 -0.079 


Total 32.0 32.5 0.896 
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Table 3: Values of % 2 , and corresponding percentage C.L. (d.o.f. = 2) of the atmospheric 
neutrino anomaly hypothesis, for the experiments analyzed in Sec. 2. Flux errors as large 
as 30%, 20% and 10% are considered. 





S flux 


= 30% 


$ flux 


= 20% 


$ flux ' 


= 10% 


Experiment 


x 2 


C.L. 


x 2 


C.L. 


x 2 


C.L. 


Kamioka sub-GeV 


12.7 


99.8 


13.5 


99.9 


15.7 


99.9 


Kamioka multi-GeV (FC) 


7.08 


97.1 


7.11 


97.1 


7.17 


97.2 


Kamioka multi-GeV (FC+PC) 


8.79 


98.8 


8.81 


98.8 


8.86 


98.8 


1MB 


6.13 


95.3 


6.53 


96.2 


14.3 


99.9 


Frejus (FC) 


0.79 


32.6 


1.26 


46.7 


2.45 


70.6 


Frejus (ALL) 


0.32 


14.8 


0.54 


23.6 


1.13 


43.2 


NUSEX 


0.14 


6.76 


0.25 


11.8 


0.48 


21.3 


Soudan 2 (preliminary) 


2.04 


63.9 


2.06 


64.3 


2.11 


65.2 
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